Anharmonic oscillator radiation process in a large cavity 
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We consider a particle represented by an anharmonic oscillator, coupled to an environment (a field) 
modeled by an ensemble of anharmonic oscillators, the whole system being confined in a cavity of 
diameter L. Up to the first perturbative order in the quartic interaction (interaction parameter A), 
we use the formalism of dressed states introduced in previous publications, to obtain for a large 
cavity explicit A-dependent formulas for the particle radiation process. These formulas are obtained 
in terms of the corresponding exact expressions for the linear case. We conclude for the enhancement 
of the particle decay induced by the quartic interaction. 
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I. INTRODUCTION 



In previous publications [1-3], a non perturbative ap- 
proach {dressed states) has been used to study systems 
that can be described by a Hamiltonian of the form, 
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(1.1) 



where the limit N — > oo is understood, the subscript 
refers to a particle approximated by a harmonic oscilla- 
tor having bare frequency w , and k = 1,2, ...N refer to 
the harmonic field modes. A Hamiltonian of the type of 
Eq.(l.l), can be viewed as a linear coupling of an atom 
with the scalar potential, and has been investigated in 
[1,2] to define rigorous dressed states which allow a non- 
perturbative approach to the time evolution of atomic 
states. An Hamiltonian of the type (1.1) has also been 
employed in [4] to study the quantum Brownian motion 
of a particle with the path-integral formalism. In the 
limit N — > oo we recover, within the framework of the 
harmonic approximation, the situation of the linear cou- 
pling of an atom with the scalar potential, or the coupling 
of a Brownian particle with its environment, after redef- 
inition of divergent quantities. In the case of the cou- 
pled atom field system, the above mentioned formalism 
of dressed states recovers the experimental observation 
that excited states of atoms in sufficiently small cavities 
are stable. It allows to give formulas for the probability 
of an atom to remain excited for an infinitely long time, 
provided it is placed in a cavity of appropriate size [2]. 

In this note we intend to generalize, making some ap- 
propriate approximations, to an anharmonic oscillator, 
the linear coupling to an environnement as it has been 
considered in the above mentioned works. In this case, 
the whole system is described by the Hamiltonian, 

N 

Hi = H(p Q , q , {p k ,q k }) + ^ KT' l ^ pa q fl q v q p q <7 , (1.2) 
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where H (po, <?0j {Pk, qk}) is the bilinear Hamiltonian in 
Eq. (1.1) and T$ pa are some coefficients that will be de- 
fined below. In Eq. (1.2) summation over repeated greek 
labels is understood. We emphasize that our problem 
is different from the situation treated in the pioneering 
papers of Refs [5,6]. We do not intend to go to higher- 
orders in the perturbative series for the energy eigen- 
statates, we will remain at a first order correction in A r , 
and we will try to see what are the effects of the anhar- 
monicity term, given by the last term of Eq. (1.2), on 
our previous results for the linear coupling with an en- 
vironnement. We notice that the anharmonicity term in 
Eq. (1.2) involves, independent quartic terms of the type 
q*, n = 0, {i} (self-coupling of the bare oscillator and of 
the field modes), quartic terms coupling the oscillator 
to the field modes and also the terms coupling the field 
modes among themselves. These terms are of the type 

qoqf, qfah ilni and of tne typ e QiiJ 1 ' with n + m = 4 ' 

for the coupling between the field modes. We intend to 
start from the exact solutions we have found in the linear 
case, and investigate at first order, how the quartic in- 
teraction characteristic of the anharmonicity changes the 
decay probabilities obtained in the linear case. 

In the section II we review some results obtained in 
the previous publications mentioned above. In section 
III we show how, making appropriate approximations, 
the dressed states approach introduced for ohmic systems 
can be used to generalize some results to an anharmonic 
oscillator coupled to an ohmic environnmcnt. 



II. THE HARMONIC SYSTEM 

The bilinear Hamiltonian (1.1) can be turned to prin- 
cipal axis by means of a point transformation, q^ — 
t^Qr, P» = t^P T ; ii = (0,{fc}), k = 1,2, ...,7V; r = 
0, ...N, performed by an orthonormal matrix T — (tV). 
The subscript \i = refers to the atom (or the Brownian 
particle) and ji = k ,k = 1,2,3... refer to the harmonic 
modes of the field (or the thermal bath) The subscripts 



1 



r refer to the normal modes. In terms of normal mo- 
menta and coordinates, the transformed Hamiltonian in 
principal axis reads, 
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(2.1) 



where the f2 r 's are the normal frequencies corresponding 
to the possible collective stable oscillation modes of the 
coupled system. The matrix elements t r are given by [1] 
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with the condition, 
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(2.3) 



We take = r)(wk) n - In this case the environment 
is classified according ton > 1, n = 1, or n < 1, 
respectively as supraohmic, ohmic or subohmic. For a 
subohmic environment the sum in Eq.(2.3) is convergent 
and the frequency uj is well defined. For ohmic and 
supraohmic environments the sum in the right hand side 
of Eq.(2.3) diverges what makes the equation meaningless 
as it stands, a renormalization procedure being needed. 
We restrict ourselves to ohmic systems. In this case, 
using the method described in [1] we can define a renor- 
malized frequency u), by means of a counterterm Suj 2 , 



Slo 2 



Slo 2 = Nn 2 



in terms of which Eq.(2.3) becomes, 
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(2.4) 



(2.5) 



We see that in the limit N — > oo the above procedure is 
exactly the analogous of naive mass renormalization in 
Quantum Field Theory: the addition of a counterterm 
— Suj 2 qQ allows to compensate the infinuty of loq in such 
a way as to leave a finite, physically meaninful rcnormal- 
izcd frequency Co. This simple renormalization scheme 
has been originally introduced in Ref. [7]. 

To proceed, we take the constant n as n = ^/2gAuj, 
A10 being the interval between two neighbouring bath 
frequencies (supposed uniform) and where g is some con- 
stant (with dimension of frequency). We restrict ou- 
selves to the physical situations in which the whole sys- 
tem is confined to a cavity of diameter L, in which case 
the environment (field) frequencies cok can be writen in 
the form 



uJk = 2kn/L, k = 1,2, ... . 
Then using the formula, 



(2.6) 
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(2.7) 



and restricting ourselves to an ohmic environment, 
Eq.(2.5) can be written in closed form, 



.LSI. ft c . u 2 L, 



(2.8) 



The solutions of Eq.(2.8) with respect to SI give the spec- 
trum of cigcnfrcqucncies Sl r corresponding to the collec- 
tive normal modes. The transformation matrix elements 
turning the system to principal axis are obtained in terms 
of the physically meaningful quantities Sl r , u>, after some 
rather long but straightforward manipulations analogous 
as it has been done in [1]. They read, 



rjSl r 



(Sl 2 r - Co 2 ) 2 + ^(3f>2 - Co 2 ) + n 2 g 2 S} 2 



L0 2 -Sll 
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(2.9) 



To study the time evolution of the system, we start 
from the eigenstates of our system, | no, ni, 712...), repre- 
sented by the normalized eigenfunctions in terms of the 
normal coordinates {Q r }, 



.(Q,*)=n 




r — i > risflst 
e ^ , 



(2.10) 



where H ns stands for the n s -th Hermite polynomial and 
To is the normalized vacuum eigenfunction. We intro- 
duce dressed coordinates q' and {g-} for, respectively the 
dressed atom and the dressed modes of the field, defined 
by [l], 



(2.11) 



valid for arbitrary L and where lo^ = { }. In terms 
of the bare coordinates the dressed coordinates are ex- 
pressed as, 




where 



a 



(2.12) 



(2.13) 



In terms of the dressed coordinates, we define for a 
fixed instant dressed states, | n n , ki, k 2 ..-) by means of 
the complete orthonormal set of functions [1] , 



.(</)= n 




(2.14) 
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where q'^ = q' , q[, Cu^ = {a), LOi}. Note that the ground 
state To in the above equation is the same as in Eq.(2.10). 
The invariance of the ground state is due to our definition 
of dressed coordinates given by Eq. (2.11). Each function 
ip Kom — W) describes a state in which the dressed oscilla- 
tor q' is in its — th excited state. Let us consider the 
dressed state | 0, 0, ...l(^), 0...), represented by the wave- 
function V'oo.. .1(^)0. ..(<?')■ ^ describes the configuration in 
which only the dressed oscillator q 1 ^ is in the first excited 
level. Then it is shown in [1] the following expression 
for the time evolution of the first-level excited dressed 
oscillator q'^, 

|o,o,...i(/i),o...) (t) = ^ /""(*)! o,o,... 1(1/), o...) (o), 

(2.15) 



where 



-iQ s t 



From Eq.(2.15) we see that the initially excited dressed 
oscillator naturally distributes its energy among itself 
and all other dressed oscillators (the atom and the en- 
vironment) as time goes on, with probability amplitudes 
given by the quantities f^ u (t) in Eq. (2.16). For /j, = in 
Eq. (2.15) the coefficients f 0v {t) have a simple interpre- 
tation: ,f ao (t) and f 0l (t) are respectively the probability 
amplitudes that at time t the dressed particle still be ex- 
cited or have radiated a quantum of frequency TiiVi. We 
see that this formalism allows a quite natural description 
of the radiation process as a simple exact time evolu- 
tion of the system. In the case of a very large cavity 
(free space) our method reproduces for weak coupling 
the well-known perturbative results [1,2]. 



III. AN ANHARMONIC OSCILLATOR IN A 
LARGE CAVITY 

The introduction of the quartic interaction term in 
Eq.(1.2) changes the Hamiltonian in principal axis from 
Eq.(2.1) into 
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^ ' ^rTpyppt ^t u 2 t ptJ'Q ri Q r2 Q r3 Q r/L , (3.1) 
r=0 

where summation over the repeated indices n , r 2 , r 3 , r 4 
is understood. In order to have a specific quartic inter- 
action, we make a choice for the coefficients T^) pa in the 



above equation, 



q-ij) _ .r.r.r.r 
1 iivpcr — h ii l v l p l a i 



(3.2) 



which replaced in Eq. (3.1) and using the orthonormality 
of the matrix V gives the Hamiltonian in principal axis, 

1 N N 
Hl = f £(i? + Q rQ 2 r) + £ KQ} ■ (3.3) 
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Since the quartic interaction, as given by Eq. (3.2), 
decouples the normal coordinates, the renormalization 
procedure will remains the same as in the absence of 
the quartic interaction. This means that the dressed fre- 
quency u) is still given by Eq. (2.4). 

Performing a perturbative calculation in A r , we can 
obtain the first order correction to the energy of the sys- 
tem, J2 r h^r, in such a way that the A r -corrected energy 
can be written in the form, 



£({A r }) = 7i^(Q r + A r e«) 



(2.16) where the first order correction is given by 
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(3.4) 



(3.5) 



For sufficiently small A r we will have quasi-harmonic 
normal collective modes having frequencies fi r + X r er . 
Accordingly we can describe approximately the system 
in terms of modified harmonic eigenstates, which can be 
written as a generalization of the exact eigenstates (2.10), 
replacing the eigenfrequencies O r by the A-corrected val- 
ues fl r + \. r e^r\ 



.(Q,t;{\ r }) = ll 




IV-^^+^K (3.6) 



From the modified harmonic eigenstates (3.6), we can fol- 
low analogous steps as in the harmonic case [1,2] to study 
the A-corrected evolution of a dressed particle, generaliz- 
ing Eq. (2.15), 

| l,0,0,...,0...)(t;{A 8 }) = £/ "(t;{A a }) 

x | 0,0,...l(i/),0...)(0) , (3.7) 

where 

f°"(t; {K}) = £ t5t-e- < ( n - +A - e - 1> )' (3.8) 

s 

For a very large cavity, from Eq. (2.9) we obtain for L 
arbitrarily large, 



to — * LirriL- 



2gQ v / 2TTc/L 



^/(n 2 - uj 2 ) 2 + ir 2 g 2 n 2 ' 



(3.9) 



from which using the definition of the coefficients 
f 0u (t;{X s }) from Eq. (3.8), and the fact that for very 
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large values of L, 2irc/L = 8u> = SCI, we have an ex- 
pression for the {A r }-corrected probability amplitude for 
the particle be still excited after an ellapsed time t, the 
quantity, 



(n 2 -uj 2 ) 2 + TT 2 g 2 n 2 



(3.10) 



From dimensional arguments, we can choose An = Af2 3 , 
where A is a dimensionless small fixed constant. With 
this choice, after expanding inpowers of A the exponential 
in Eq. (3.10), we obtain to first order in A the amplitude, 



15Af d 



(3.11) 



where f 00 (t) is the probability amplitude for the har- 
monic problem, that the particle be still excited after a 
time t [1,2]. For ^ < Q (situation that includes weak 
coupling, g <C u)) and for a very large cavity, f 00 (t) is 
given by [3] , 

f°(t) = (l- ^ e — *-9*/2 + l J(t)j (3 12) 



probability (3.14), for A = 1/50. We see clearly the en- 
hancement of the particle decay induced by the quartic 
interaction. 



D.4-, 



P[t] 0.2 




FIG. 1. Plot on the vertical axis, commonly named P[t], of 
the A-corrected probability (3.14) (dashed line) and the har- 
monic probability Eq.(3.15) (solid line), for u = 4.0x 10 14 /s 
and g = auj. a is the fine structure constant, a — 1/137 and 
time in units of 10 _13 s 



where, 



J(t) = 2g 



dy 



y 2 e-y t 



(y 2 + u) 2 ) 2 — Tr 2 g 2 y 2 



i t 3 
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From Eq. (3.11) we can obtain at order A, the proba- 
bility that the particle remains in the first excited state 
at time t, 
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l/ 00 (t;A)| 2 = |/ 00 W| 2 + ^||/ 00 W| 2 



(3.14) 



We know that |/ 00 (i)| 2 is a decreasing function oft, what 
means that the derivative of this function with respect to 
t is negative. Therefore, since A > 0, we conclude from 
Eq. (3.14) that |/ 00 (t;A)| 2 is smaller than the harmonic 
probability |/ 00 (t)| 2 . Indeed we know from Refs. [2,3] 
that for large t (t >> 4) we have, 



irw = (i+ 
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4k 2 
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(3.15) 



which is a rapidly decreasing function of t. 

In Fig.l we plot on the same scale the A-corrected 
probability (3.14) and the harmonic probability in 



Eq.(3.15), for ui = 4.0 x 10 14 /s and g 



where a 



is the fine structure constant, a = 1/137 and time is 
rescaled as t x 10~ 13 s. The solid line is the harmonic 
probability (3.15) and the dashed line is the A-corrected 
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